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Abstract 

A geometric framework for describing and solving time-dependent implicit differential equa- 
tions F(t, x, x') = is studied, paying special attention to the linearly singular case, where 
F is affine in the velocities: A(t,x)x' = b(t,x). This framework is based on the jet bundle 
of a time-dependent configuration space, and is an extension of the geometric framework of 
the autonomous case. When A is a singular matrix, the solutions can be obtained by means 
of constraint algorithms, either directly or through an equivalent autonomous system that 
can be constructed using the vector hull functor of affine spaces. As applications, we con- 
sider the jet bundle description of time-dependent lagrangian systems and the Skinner-Rusk 
formulation of time-dependent mechanics. 
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1 Introduction 



It is known that the dynamics of an autonomous mechanical system given by a lagrangian 
function L: TQ — > R is described by equations of motion of the type W(q, q) q = a(q, q); writing 
them as a first-order differential equation on TQ, we obtain 

A{x)x = b(x), 

where the matrix A is singular if the lagrangian is not regular (see for instance |Car 90| ^ . Indeed, 
such differential equations appear in many other applications as control theory, circuit theory, 
engineering, . . . — see some references in GMR 04 . 

A general framework to deal with differential equations of this type was given in |GP 91 j and 
developed in |QP 921 FGP 02) : the so-called linearly singular systems. Since the matrix A may be 
singular, the system may not have solutions passing through each point of TQ, and the solutions 
may not be unique. To solve this problem, a constraint algorithm has to be performed. If the 
algorithm ends, one obtains a submanifold of TQ where there exist solutions of the system. 
In fact, this algorithm is a generalization of the presymplectic constraint algorithm GNH 78 , 
which is in its turn a geometric version of the Dirac-Bergmann theory of constrained systems. 

The main purpose of this paper is to extend these results to the time-dependent case. We 
study the geometric framework of time-dependent first-order implicit differential equations, 

F(t,x,x) = 0, 

and the linearly singular case, when F is affine in the velocities, 

A(t, x)x = b(t, x), 

where A is a possibly singular matrix. 

Our model for the time-dependent configuration space, rather than a trivial product M = 
R x Q, is a fibre bundle Sau89 p: M — > R, where the base R contains the time variable. 
Such an M is isomorphic to a product R x Q, but in practical applications there may not be 
a privileged trivialization, and a possible extension to deal with field theory of course should 
not be based on a trivial bundle. Some references about time-dependent lagrangian systems are 
|EMK 911 IHF931 ICLM 9jj in the product case and |Kru 971 IMPL031 IMS 981 ILMMMK 02j in the 
fibre bundle case; see also references therein. Time-dependent systems in general are studied in 
many books, as for instance |AM 781 lOlv 93] . 

The main difference between the autonomous and the non-autonomous case is the usage of 
tangent bundles and jet bundles respectively. To describe an autonomous differential equation 
on the configuration space Q we use the tangent bundle TQ, which is a vector bundle. To 
describe a non-autonomous differential equation on the time-dependent configuration space M, 
we use its jet bundle J 1 /), which is an affine bundle over M. To describe a linearly singular 
equation on M we will use affine morphisms defined on this affine bundle. 

We will also propose a constraint algorithm for time-dependent singular systems. This 
algorithm is the natural generalization of the algorithm for the autonomous case — both in 
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the general implicit RR94, M MT95] and linearly singular cases |GP911 IUP92] . The case of 
an implicit equation in a product M = R x Q has already been discussed in |Del 04| . It is 
worth noting that constraint algorithms for some time-dependent systems have been described 
in several papers, as for instance |CF 931 IHLM9! I1LMM 9~9l ILMM 961 ILMMMR 02l | Vig 00| . All 



these systems are of linearly singular type, so they are included within our framework. 

When studying a time-dependent differential equation, sometimes it is useful to convert it 
into an equivalent time- independent one. This is even more interesting for implicit equations; 
for instance, the constraint algorithm for the autonomous case is easier to implement than for 
the non-autonomous case, because of the fact that vector fields instead of jet fields are used to 
obtain the constraint functions. 

Therefore, we will examine the possibility of associating an autonomous linearly singular sys- 
tem with a time-dependent one, so that the solutions of both systems will be in correspondence. 
Essentially, we use the canonical inclusion of J 1 /? into TM. In order to perform this association, 
we propose two different strategies. One possibility is to choose a connection on the jet bundle 
to induce a splitting of the tangent bundle. The other possibility, which does not make use of 
any choice, is based on the notion of vector hull of an affine space or an affine bundle. The main 
idea is that any affine space A can be canonically embedded as a hyperplane in a vector space A 
— the vector hull of A; with this immersion affine maps can be homogenized, that is, converted 
into linear maps. 

As we have pointed out at the beginning, our main motivation comes from Euler-Lagrange 
equations and mechanical systems, where equations of motion are of second order. So it is 
interesting to extend the preceding study to second-order implicit and linearly singular equations: 

F(t, x, x, x) = 0, A(t, x, x)x = b(t, x, x). 

As applications of the formalism, we give two descriptions of time-dependent mechanical 
systems, in the form of time-dependent singular lagrangian systems and in the mixed velocity- 
momentum description (sometimes called Skinner-Rusk formulation [Ski 831 SR83 ) of time- 
dependent mechanics DMH_Q2j- As a concrete example, we also study a pendulum of variable 
length. 

The paper is organized as follows. In section 2 we study the geometric formulation of 
time-dependent differential equations, either in the implicit and in the linearly singular case. 
In section 3 we describe constraint algorithms for both cases. In section 4 we present two 
constructions of an autonomous system associated with a given time-dependent system. The 
extension to second-order equations is presented in section 5. Applications to singular lagrangian 
mechanics are presented in section 6, and section 7 is devoted to an example. Finally, there is 
an appendix about vector hulls of affine spaces and affine bundles. 



2 Time-dependent systems 

In this section we discuss first-order time-dependent singular differential equations. As a model 
of time-dependent configuration space, we take a fibre bundle p: M — > R over the real line 
(though more general settings could be also considered). 
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The appropriate geometric framework to deal with derivatives is that of jet bundles, so we 
will begin by giving some facts and notation about them — see for instance |Sau 89j . 



2.1 First-order jet bundles 

We denote by J 1 /? the first-order jet manifold of p. Its elements, called 1-jets, are equivalence 
classes of local sections of p: two sections are equivalent at a point when they are tangent. We 
denote by j}£ the 1-jet of a section £ at t. 

J 1 /? is a fibre bundle over M and over R, with canonical projections shown in the diagram: 




R 



We denote by t the canonical coordinate of R. If (t, q % ) are fibred coordinates on M, then J 1 /? 
has induced coordinates (t,q l ,v l ). 

The bundle p± : o: J 1 /? — > M is an affine bundle modelled on the vertical bundle of p, Yp — > M. 
Recall that the vertical bundle Yp is a subbundle of TM, given by Yp = KerTp. The elements 
of Yp are the tangent vectors of M which are tangent to the fibres; equivalently, Yp = {v € 
TM | i v dt = 0}. 

There is a canonical embedding r.^p — > TM, defined as i(j|£) = In local coordinates, 
t{t, g*', «*) = (t, g 4 , 1, u 4 ). Notice that t(J x p) = {d£ TM | z„dt = 1}. 

2.2 Implicit systems 

In general, a (time-dependent) implicit differential equation is defined by a submanifold X> C J 1 p. 
A local section £: / — * M of p is a solution of the differential equation if 

j 1 ^) e f (i) 

for each i. If the subset P is locally described in coordinates by some equations F a (t, q l ,v l ) = 0, 
then the differential equation reads F a (t, C(t)) = 0. 

Suppose that T> is the image of a jet field, that is, of a section X: M — ► J 1 /}. Then the 
solutions of the differential equation are the integral sections of X, which are the solutions of 
the explicit differential equation 

3 1 Z = Xo£. 

In coordinates, if X(t,q l ) = (t,q t ,X^(t,q' t )), the differential equation reads = X l (t,^(t)). 

Consider again T> C J x p. Given a jet field X, its integral sections are solutions of the implicit 
equation defined by T> iff 

X(M) C V. (2) 

So, in a certain sense, solving this equation is equivalent to solving the implicit equation l(T|). 

For an explicit differential equation there always exist solutions, and each initial condition 
x € M defines a unique maximal solution. For an implicit differential equation existence and 
uniqueness may fail; in this case one is lead to study the subset of points covered by solutions, 
and the multiplicity of solutions. 
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2.3 Linearly singular systems 

A (time-dependent) linearly singular system on M is defined by a vector bundle n: E — > M and 
an affine bundle morphism A: J 1 /? — ► E: 

JV — £ (3) 

Pl,0 

M 

For the sake of brevity, we will refer to this linearly singular system simply as A. 
The system A has an associated implicit system given by 

P = ^ 1 (0)CJV (4) 

A local section £:I — > M is a solution of D, equation (J2), iff it is a solution of the linearly 
singular differential equation 

Ao^ = 0. (5) 
In local coordinates, the bundle morphisms are given by 

7r(t,q\u a ) = (t,<f), A(t, = (t^^AJ^q'y +c a (t,q% 

thus the differential equation reads 

A°(t,e(t))e(t)+ c <*(t,e(t)) = o. 

As before, it may be convenient to describe the solutions of the differential equation as 
integral sections of jet fields. A jet field X: M — > J 1 /? is a solution jet field of V, equation ©, iff 

AoX = 0. (6) 

Then its integral sections are solutions of the differential equation defined by A. 
Locally, X(t,q l ) = (i, q 1 , X l (t, q 1 )) is a solution jet field of A when 

A^(t,q i )X^t,q i ) + C a (t,q i )=0. 

Let us remark that, instead of a vector bundle, we could have considered an affine bundle 
E with a section b and an affine bundle morphism A: J 1 /? — > E. This slight generalization does 
not seem too relevant for applications, and indeed the section b endows E with a vector bundle 
structure. 

3 Constraint algorithm 

In general, an implicit system does not have solution jet fields, and does not have solution 
sections passing through every point in M. We want to find a maximal subbundle p': M' — > R 
of p (over R for simplicity, but more general situations could occur) where there exist solution 
jet fields X: M' — > J 1 ^' and solution sections £: I — ► M' through every point in M' . 

To this end, we can adapt the constraint algorithms of the time-independent case, both for 
implicit systems (RR 94j MMT95] and linearly singular systems |GP 911 RlP 92j . to the time- 
dependent case. A constraint algorithm for a time-dependent implicit equation in a product 
M = RxQ has been recently discussed in |Del 04| . 
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3.1 Implicit systems 

Let T> C J 1 /? be an implicit system. We say that a 1-jet y € T> is integrable (or locally solvable) 
if there exists a solution £: / ^ M of 2? such that j 1 ^ passes through y. One of the purposes of 
the constraint algorithm is to find the set T>i n t of all integrable 1-jets. 

If a solution passes through a point x £ M, necessarily x belongs to the subset 

M (1) :=p lfi (V). (7) 

Denote by /jmiMm — > R the restriction of p to Mm. To proceed with the algorithm, we will 
assume that pm is a subbundle of p. In this case, the inclusion imiMm M has a 1-jet 
prolongation, j^miJ Vi) <— > J 1 /). By means of this inclusion, we can define 

£>(l) := JV(l)HP. (8) 

Since the solutions of P lay on Mm , the integrable jets of P must be contained in . If this 
is a submanifold, we have obtained a new implicit system, now on Mm. 

This procedure can be iterated: from Mm = M and X>(q) = £>, and assuming that at 
each step one obtains subbundles and submanifolds, one may define Mo\ := pi.o(P(j-i)) an d 
T>u\ := J 1 p(j) fl T>u-yy The algorithm finishes when, for some k, we have M/ k+1 \ = M/ k y In 
this case, since pi,o(P(fc)) = M( fc ), ^ we suppose for instance that the projection — > Mo.\ is 
a submersion, we have that Pi nt = 

3.2 Linearly singular systems 

Let A: J 1 /? — ► £7 be a time-dependent linearly singular system as described in section 2. We 
can proceed by applying the preceding algorithm for implicit systems, and also by adapting the 
algorithm for time-independent linearly singular systems. 

So we begin with V = .4 -1 (0). As before, the configuration space must be restricted to 
Mm := pifl(T>), which can also be described as 

M (1) = {x £ M | X e 1m Ax}; 

note that X £ Im^ is the necessary consistency condition for @ to hold on a given point 
x 6 M. 

As above, we assume that pm: Mnj —* R is a subbundle of M. 

Let us restrict all the data to Mm: .4m := A\ji p , Sm := -^Im (1) ' anc ^ ^(i) := ■ ^° 

we obtain a linearly singular system on Mm: 

J V(i) — ^(D (9) 

(P(l))l,0 

M (1) 

It is clear that the implicit system defined by Am coincides with the implicit system 2?m 
obtained above, (JSJ), that is, 

A^ ) (o) = J 1 p {1) nv. (10) 
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Thus, if we assume for instance that each is a subbundle of M^_^ , we obtain a constraint 
algorithm for the linearly singular case: 

:= {x € M (i _ 1} | 0^ e Im^-i))*}, 



JV(i) 



(0 



M, 



(<) 



7T(j) := TT(i-l) 



(i) 



When the algorithm finishes, we arrive to a final system which is integrable everywhere. 




4 From non-autonomous to autonomous systems 

It is usual to convert a time-dependent system into a time-independent one by considering the 
evolution parameter as an new dependent variable. From a geometric viewpoint, this is easily 
done with an implicit system D C J 1 /} by means of the canonical inclusion i: J 1 /? — > TM: its 
image E = i(D) is a submanifold of TM, so it defines an autonomous implicit equation. The 
equivalence between both equations is immediate: 

Proposition 1 A map £: R — > M is a solution section of the time- dependent system D iff it is 
a solution path of the autonomous system E such that p(£(io)) = *o f or an U arbitrarily given t$M 



Now let us focus on the case of a time-dependent linearly singular system A: 



Pl,0 




(11) 



We will also relate this system to an autonomous one. The main motivation for finding such a 
relation is that the constraint algorithm described in the preceding section is easier to implement 
in the autonomous case. The reason is that, instead of jet fields, vector fields can be used to 
obtain constraint functions that define the submanifolds in the constraint algorithm, as will be 
shown later on. 

Two constructions to achieve our goal will be proposed. In the first one, we use a connection 
to define a complement of V p in TM. In the second construction, we use the vector hull functor 
described in the appendix to define vector bundles and morphisms from affine bundles and 
morphisms. 
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Previously, we shall recall some definitions concerning the autonomous case |GP 92j . An 
autonomous linearly singular system on a manifold A is defined by a vector bundle it: F — > A, 
a vector bundle morphism A: TN — > i* 1 , and a section 6: A — > i* 1 : 

A 

Taking local coordinates i; 2 ) on TA and (a; 1 , u Q ) on F, the local expressions of the maps are 

Tr(x i ,u a ) = {x i ), b(x i ) = {x i ,b a (x i )), A(x\x i ) = {x i ,A?{x i )x j ). 
We say that a path 7: / — > TV" is a solution path if 

A o 7 = b o 7. 

Locally, (7(^)7-^ (t) = b a {^f{t)). A vector field A E X(A) is a solution vector field when 

ioI = i); 

in coordinates, Aj(x)X J (x) = 

Let us roughly recall — see GP 92 for a detailed description — how the constraint algorithm 
is explicitly carried out in this autonomous case, that is, how the consecutive constraint manifolds 
are effectively computed. First of all, it can be seen that the primary constraint submanifold 
M\ = {x E M I b(x) E ImAj} is locally described by the vanishing of the functions <p a := {s a , b), 
where (s a ) is a local frame for Ker'A C F* . The constancy of the rank of A is needed here to 
ensure that this procedure works. 

We have that a vector field X in M, in order to be a solution of the system, must satisfy the 
equation Ao X ~ b. Vector fields satisfying this condition always exist and are called primary 

Mi ' 

vector fields. Given one primary vector field Ao, the others have the form X ~ X$ + J2/i 

where jf* are functions uniquely determined on M\ and (r«) is a local frame for Ker A. 

A primary vector field X can be a solution of the system only if it is tangent to M±, so 
we obtain the equation, for every constraint (f> a , (X ■ (j) a ){x) = 0, for x E Mi, or, equivalently, 
(Ao- 4> a )(x) + J2 u(r 'ft' 4> a ){x) j '(x) = 0. These equations may provide new constraints that define 
the secondary constraint submanifold M2, and may also determine some of the functions / M . 

This procedure can be iterated until we determine which primary vector fields are solutions 
of the system, and we obtain the final submanifold where they are defined. 

4.1 Jet field construction 

Consider the linearly singular system given by Let us choose an arbitrary jet field T: M — > 
J x p. This jet field T induces |Sau 89 in a natural way a connection T on the bundle p and a 
splitting TM = Vp © Hr of the tangent bundle, with projections vr and hr- The coordinate 
expressions are: 

T(t,q i ) = (t,q i ,r(t,q i )), 
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Vr(t, g*; i, <f ) = (t, g*; 0, g 4 - f g*)), fc r (t, g'; i, <f ) = (t, g*; t, f P (f, g*)). 
From this we can define a section of tt, 

6 r := -AoT:M -> E, 

and a vector bundle morphism, 

A T := Aov r :TM -> E 1 , 

where *4: Vp — > £7 is the vector bundle morphism associated with the afhne map A. 
With these objects we can construct a time-independent linearly singular system: 

TM — £0R (12) 

TM / 

M 

This system is equivalent to the time-dependent system in the sense of the following 
proposition. 

Proposition 2 Consider the time- dependent system given by Ml\) . Given any jet field V: M — > 
J 1 /?, we have: 

i) A map £: R — > M is a solution section of the time- dependent system 111]) if, and only if, it 
is a solution path of the autonomous system M/Jjl such that p(£(to)) = ^0 f° r an D arbitrarily 
given to. 

ii) A map X: M ^ J 1 /} is a solution jet field of the time- dependent system h 11]) if, and 
only if, considered as a vector field in M , it is a solution vector field of the autonomous 
system \lty) . 

(Note that we are using the embedding i: J 1 /) — > TM defined in section 2 to identify jet fields as 
vector fields.) 

Proof. It is immediate, taking into account the local expressions of the equations defined re- 
spectively by both systems: 

• Af(t,^)vj +c a (t,q i ) = 



AfM)($-iTi(t,tf)) = -Af(t,q i )TS(t,q i )-d>(t,q i ) 
i = 1 



9 



4.2 Vector hull construction 



Here we will apply the vector hull functor described in the appendix. The affine bundle morphism 
A in induces a vector bundle morphism A between the vector hulls of J l p and E: 



A 
E(- 



■Jip 
A 
~E 



The section of ir: E — ► M also induces a section of n : E — > M, defined by = i o 0. Recall 
that with the identification E = E x R, we have = (0, 1). 

Using the canonical identification of J x p with TM, we can construct the following linearly 
singular system: 



TM 

M 



A 



E 



(13) 




This system is equivalent to the time-dependent system IjllJI in the sense of the following 
proposition. 

Proposition 3 Consider the time- dependent system given by Ul\) . 

i) A map £: R — > M is a solution section of the time- dependent system hll\) if, and only if, 
it is a solution path of the associated autonomous system such that p(£(to)) = to for 
any arbitrarily given to. 

ii) A map X:M^> J 1 /? is a solution jet field of the time- dependent system Ml\) if, and only if, 
considered as a vector field in M, it is a solution vector field of the associated autonomous 
system kltft) . 

Proof. Again we can prove the result in local coordinates, where the equations of the systems 
(fTT|) and (|T5|) read, respectively, 

• c a (t,q i )+A^(t,q i y=0 
f c a (t,q i )i + Af(t,q i )q^ =0 

We have used the coordinates (u°,u a ) on E induced by the affine frame (0;u a ) on E. ■ 



5 The second-order case 

The preceding results can be extended to consider higher-order implicit and linearly singular 
differential equations. Of course, the most important case, due to its applications to mechanics, 
is that of second-order equations, to which we devote this section. Before proceeding, we need 
some additional facts about the second-order jet bundle of a fibre bundle p: M — > R. 
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5.1 Second-order jet bundles 

The jet space 3 2 p is a fibre bundle over J 1 /?, M and R. The canonical projections are: 

M 
p 
R 

The bundle p2 1: J 2 p — > J 1 /? is an affine bundle modelled on the vertical bundle of pi o, Vpi o — > 

There exists a natural embedding j: J 2 p — > J 1 /?!, defined by j(j 2 £) = jj (j 1 ^)- I n addition to j, 
we can use the embedding t: J 1 pi — > TJ 1 /? to include the second-order jet bundle into a tangent 
bundle: 

j 2 p jVi tjV 

Consider natural coordinates (t,q i ,v i ,a i ) on J 2 p. With them, the local expression of this 
composed embedding k = i o j reads 

/e(t, a*) = 1,^,0*). 

This shows that the image of J 2 p by the embedding can be expressed as 

k(J 2 P ) = {w G TJ 1 p | i^dt = 1, S^iu) = 0}. 

Here there appears another relevant operator, the canonical vertical endomorphism S of TJ 1 /?, 
whose local expression is 

S= (dg* - «Mt) ® A ; 

note also that Im(5) = Vpi o- 

The Cartan distribution of J 1 ^ is the kernel of the vertical endomorphism S of T(J 1 p). We 
denote it by C7pi,o- Locally, we can describe Cp^o as the distribution generated by the n + 1 
vector fields { + f*^r, ^r}- We have an exact sequence 

> Cpi fi T(JV) Vpi, (14) 

5.2 Second-order implicit and linearly singular systems 

Similar to the first-order case, a second-order implicit differential equation is defined by a sub- 
manifold T> C J 2 p. A local section £: / — > M is a solution of the differential equation if j 2 £(t) G X> 
for each i. In coordinates, this equation can be expressed as F a (t, = 0- Like in 

the first-order case, if V is the image of a section X: J 1 /? — > J 2 p, the equation can be written in 
normal form. 



P2,0 
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Now let us consider the linearly singular case. A time- dependent second-order linearly singu- 
lar system is denned by a vector bundle tt: E — > 3 1 p and an affine bundle morphism A: J 2 p — > -E: 

JV — £ (15) 




P2,l 

Its solution sections are sections £ of p such that 

-4°j 2 £ = 0. 

Locally this reads 

A«(t,e(t),e(t))e(t) + c a (t,e(t),e(t)) = o. (ie) 

A second-order jet field, that is, a section X of P(2,i) : J 2 / 9 ~ * J 1 / ) is a solution jet field if 

iol = 0; 

locally this reads Af(t,tf,v i )X'(t,q i ,v i )+(? l (t,q i ,v i ) = 0. 

In a similar way of what we did in the previous section, it is interesting to convert the 
singular system given by (|15[) into a first-order autonomous linearly singular system. As before, 
we present two constructions of this. 

5.3 Jet field construction 

Choose an arbitrary second-order jet field T: J l p — > J 2 p. Again this determines a splitting of 
the tangent bundle of J 1 /) as a direct sum TJ 1 /? = Vpip © Hp, with associated projections vr 
and hr- Now we define a section of it 

b r ■= -AoT-.Pp-tE 

and a vector bundle morphism 

A T := Aov r :T3 1 p^ E. 

With these definitions, we obtain an autonomous first-order linearly singular system on the 
manifold J 1 /?: 

TJV AreSed *> £ @ Vpi,o®R (17) 




-®0©1 

A result quite similar to Proposition |21 can be formulated, in the sense that this system is 
equivalent to the original time-dependent second-order system (|15j) . This can be readily seen by 
comparing the local expression (|16|) with that of the equation defined by system (|17jl. We skip 
the details. 
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5.4 Vector hull construction 



As opposite to the first-order case, the vector hull of J 2 p can not be identified with a tangent 
bundle, but rather with a tangent subbundle. As it is shown in the appendix, J 2 p can be 
identified with the Cartan distribution Cpio on J 1 p. Then, as in section POl if we homogenize 
the system (|15|l we obtain the following: 

Cpi,o^ E (18) 



Til 




This is an autonomous linearly singular system on J 1 ^, except for the fact that there is only 
a subbundle Cpifi C TJ x p instead of the whole tangent bundle. The interpretation of this 
system is the same as in the ordinary case, but with the additional requirement that, for a path 
77:/ —> J 1 /), its derivative r\ must lie in Cpifl — which is a natural condition if r\ has to be the 
lift j 1 ^ of a section of p. In coordinates, if rj = (t, q % ,v % ), this requirement amounts to 

c? = iv\ 

Assuming that fj C Cpifl, the local equations for the path 77 to be a solution of the system (|18j) 
are 

!. 1 

\ tc a {t,q\v l )+v^A^(t,q\v l ) = 

Comparing these three equations with equation (|TE|) . we see that systems (fT5|) and (fT%|) are 
equivalent. 

Finally, we will show that the system (|18|) can be expressed as a linearly singular system, 
provided we have an appropriate extension of A. Since E is a vector bundle, we have a canonical 
identification E = £®R, and the vector extension A can be written A = AC © dt (see the 
appendix). Suppose that we have an extension A:T{J l p) — > E of the map JC:Cp\jo — ► E\ in 
some applications (see next section) there exists a natural extension A. Then the system ()18j) 
can be described as the following linearly singular system: 

tjV Aedms > ffeR@V/>i, (19) 

jy ' 

The only thing to be noted is that C pio is the kernel of S, see (|14|) . 



6 Some applications to mechanics 
6.1 Non-autonomous lagrangian systems 

Let p: Q — > R be a fibre bundle modelling a time-dependent configuration space. A lagrangian 
system on Q is determined by a lagrangian function L: J 1 /? — > R. 
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The vertical endomorphism S of T(J 1 p) allows to construct the Poincare-Cartan forms 

@L = t S o dL + Ldt e n 1 (3 1 p), 

fl L = -dQ L e J7 2 (jV)- 

By contraction, this one defines a morphism f^iTp 1 p) — > T*(J 1 p). 

In coordinates it, q, v, a) of J 2 p, the Euler-Lagrange equation can be written dp/dt = dL/dq, 
where p = dL/dv. Now, given a vector field X G X(J 1 p), we can compute 

i x n L = ((X-q) - v(X-t)j dp - ((X-p) - ^(X-t)) dq + [{X-p)v - ^{X-q^j dt. 

d d d 

Consider the case where X is a second-order vector field, X = — + v— — h A(t, q, v) — , which 

at oq ov 

amounts also to 



i x dt = l, SoX = 0. 



(20) 



Then the preceding expression simplifies to ix^L = y{X-p) — ~~q~J (vdt—dq), and so the integral 
curves of X are a solution of the Euler-Lagrange equation iff 



dL\ 



ix^L = 0. 



(21) 



Now recall the affine inclusion k: J 2 p =— > TJ 1 /?, which identifies jet fields J 1 /) — > J 2 /0 with 
second-order vector fields on J 1 /?. From the preceding discussion, it is clear that the lagrangian 
dynamics may be described by the following second-order linearly singular system on Q: 



J 2 P 

p2,l 

jV 



S!f,OK 



T*J L p 



(22) 



3h 



Using the vector hull construction described in the preceding section, we can convert this system 
into a first-order autonomous system on J 1 /?: 



Tjy 



n L ®&t®s 



T* 3 1 p R © V/91,0 



'oeieo 



and note that its equations of motion are precisely (|2Uj) . (|21|). If the lagrangian is regular, then 
this linearly singular system is regular; otherwise, the system is singular and the constraint 
algorithm for linearly singular systems can be applied to obtain the dynamics. 

Finally, let us note that there are other equivalent descriptions of the dynamics. For instance, 
instead of Ql, we could have used the Euler-Lagrange form on J 2 p. We omit the details. 
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6.2 Skinner-Rusk formulation of time-dependent mechanics 

A mixed lagrangian-hamiltonian formulation of time-independent mechanics was studied geo- 
metrically in a series of papers by Skinner and Rusk S ki83llSR83j . Recently, the time-dependent 
case has been studied in [CMC 02] . We will show how this can be described in our formalism. 

Our starting point is a fibre bundle p: Q — > R and a lagrangian function L: J 1 /? — > R. In 
this formulation, the dynamics is represented by a first-order system on the manifold M := 
T*Q X Q J 1 / - Denote the several projections as in the following diagram: 

f >> 

M :=T*Q X Q JV-^JV 




We can define the following function on M: 

U = (pr l5 pr 2 ) - pr^L, 

where ( , } denotes the natural pairing between vectors and covectors on Q, and the 2-form on M 

fin = pr* luq — dTC A dt, 

where ojq is the canonical symplectic form on T*Q; it defines a morphism 0^: TM — >• T*M. 

With these definitions we can write the equations of the dynamics in the Skinner-Rusk 
formulation, which are 

iz^H = 
i z dt = l ' 

for a vector field Z £ X(M). These equations are equivalent to the time-dependent linearly 
singular system on M defined by the following diagram: 




7 An example: a simple pendulum of given variable length 

Consider a simple pendulum whose length is given by a time-dependent function R(t). Its 
equation of motion can be written as [HR 92 



x 

y 



x 2 + y 2 



-TX 



-TV ~ 9 
R 2 (t) 
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where g is the gravitational acceleration and rR(t) is the string tension per unit mass. 

This system can be described as a time-dependent linearly singular system in the following 
way. Take M := R 6 , with coordinates (t,x,y,v x ,v y ,T), as a configuration manifold fibred over 
R, with coordinate t. The product M x R 5 is a trivial vector bundle over M, and the affine 
bundle morphism A: J 1 /? — » M X R 5 defined by 

A(x, y, v x , v y ,f)p = (x - v x ,y - v y , v x + tx, v y + ry + g, -{x 2 + y 2 - R 2 (t))) p , 

where p = (t, x, y, v x ,v v ,t) G M, models the system. 

Choosing the connection f = dt © Jj, as described in section we can convert this system 
into an autonomous linearly singular system on M, which can be written as 



1 


1 











o \ 


/ 


i 








\ 








1 













X 






Vy 













1 










y 






— TX 
















1 







Vx 






-ry-g 
























Vy 






x 2 + y 2 -R 2 (t) 




1 1 














o / 


I 


T 


) 


V 


1 


/ 



We solve this autonomous linearly singular system by means of the constraint algorithm for 
the autonomous case that is sketched in section 4. In this case, three steps are needed to solve 
the system. We give here only the constraint functions eft 1 and manifolds Mi obtained at each 
step: 

= x 2 + y 2 -R 2 (t), 
Mx = {0 1 = 0}, 

and the possible solution vector fields are of the form 



1. 



X 



d 



— — h v x — — h v v — tx— (ry + g)— 1- f-^—, 

:i/ 1 oy ov x ov y :] " 



d d d 
Mi dt ^ dx y dy dv, 

where / G C°°(Mi) is a function to be determined 
eh 2 = X ■ 6 l ~ xv x + yv v - RR', 

Mi 

M 2 = {0 1 = ^ 2 = 0}. 

03 = x ■ cj) 2 ~ v 2 x + v 2 - tR 2 - (RR" + (R r ) 2 ), 

M 2 y 

M 3 = {(j) 1 = <p 2 = cj) 3 = 0}. 

jA _ v . i 3 



a 



d_ 



— X ■ cf) 6 ~ -4tRR' - 3v y g - (RR'" + 3R'R") - fR 2 . 

The equation (ft 4 = determines the function / G C°°(M3) 
and the system is solved. 



so it is not a new constraint 



Using polar coordinates (t, r, ip, v r , v v , r), defined by 



x = r cos (f 

y = r sin <p 

v x = v r cos ip — v v r sin ip 

Vy = v r sin ip + v v r cos ip 
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we see that the final submanifold M3 is diffeomorphic to R x TS 1 , and it is embedded in M by 
(t, tp, v v ) (t, ¥>, fl'(t), (u*fl(t) 
is 



dt 



gsincp — R" (t)) / R(t)) . The (unique) solution vector field 
d 2R'(t)v<p + g cos ip d 



R(t) 



dv 



Appendix: vector hulls 
Vector hulls of affine spaces 

Though affine geometry is well-known, it presents some less known canonical structures. The 
basic fact is that any affine space has a canonical immersion, as a hyperplane, in a vector space; 
with it, affine maps can be understood as linear maps. We will not go into the details of this 
construction, nor into its several applications — see for instance : BB 75llBer 871 IGM 04l lMMS 02 , 
instead we will describe some of the basic facts. 

If A is a real affine space, we denote by A its associated vector space. We denote by A(A, B) 
the set of affine maps between the affine spaces A and B. For any affine map /: A — > B there is 
an associated linear map / G C(A, B). 

The vector hull of A is a vector space A together with an affine map j:A — > A, and such 
that the following universal property holds: for every vector space F and affine map h: A — > F, 
there exists a unique linear map h": A — ► F such that h = h'o j. It turns out that j is an affine 
immersion, with j(A) a hyperplane in A not containing 0. 



A ' 




A 

1 

I K 
Y 

F 



Notice that such a hyperplane can be described as the set w (1), with w. A — > R a unique 
linear form. This also identifies A with w _1 (0). We can gather all this information in a diagram: 



A 



A- 



A 



R 



From the universal property, it is clear that the assignment h \— > 
A(A,F) ~ C(A,F), and in particular we have A(A,TL) ~ A*. 

Given an affine map f: A — > B, there is a unique linear map /: A 
We will call it the vector extension of / for obvious reasons. 



K is indeed an isomorphism 
-> B such that f oj A = j B °f- 



A 



f 



B 



OA 



Ob 



B 
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Now the assignment / i— > / is an affine inclusion A(A, B) B), 

Now let us describe the particular, but important, case where where the affine space is a 
vector space F. Then we have a canonical identification F = R x F, with the inclusion jp{u) = 
(1, it), and wf(X, u) = A. Given an affine map h: A — > F, its vector extension h: A —> F = Hx F 
is given by h = (wa, h"). 

More particularly, for an affine map h:E^F between vector spaces, h(u) = ho + h\ ■ u, we 
have that K: E x R — > F is given by /f (A, u) = /ioA + h\ ■ u, and the vector extension of h is 
h(X, u) = (A, hoX + hi ■ u). 

Finally, let us put coordinates everywhere. Consider a point eo = ao € A and a basis (ej)j G j 
of A Then, with the appropriate identifications, every point in A can be uniquely written as 
x = x°eo + x l ei. The point x belongs to A iff x° = 1, and belongs to A iff x° = 0. With these 
coordinates, the vector extension of an affine map y 3 = c? -\-T\x % is the linear map y u = T^x^, 
withr J = c^,T{j = l, T° = 0. 

Vector hulls of affine bundles and jet bundles 

All that we have done up to now with affine spaces can be formulated for affine bundles in an 
analogous way MMS 02 . In this case, the starting point is an affine bundle A — > M modelled on 
a vector bundle A — > M. Without going into technical details, the vector hull of A is the vector 
bundle A —* M whose fibres are the vector hulls of the fibres of A (which are affine spaces). All 
the affine and linear maps that we considered before are now affine and vector bundle morphisms 
over the identity map on M. 

For our discussion it is particularly useful the following fact: if we have an exact sequence 

of vector bundle morphisms *~ A — — ■*■ W — M x R >■ , then w _1 (l) is an affine 

bundle modelled on A (so it is isomorphic to A), and W is canonically isomorphic to the vector 
hull of w-^l). 

We will apply this to find the vector hulls of the affine bundles that play a role in this paper, 
that is, the jet bundles p\$: J 1 /? — > M and p2,i' J 2 p — > J 1 /?, when p: M — > R is a fibre bundle over 
the real line. Recall sections 2.1 and 5.1 for the definition and properties of these jet bundles. 

From section 2.1 it follows that the sequence of vector bundle morphisms 

^ Vp C ^ TM M x R ^ 

is exact, where dt denotes, by abuse of notation, the contraction of a tangent vector with the 
1-form dt. Furthermore, it is also seen in section 2.1 that J 1 ^ can be canonically embedded 
into TM, and the image of this embedding is just dt _1 (l). Therefore, the vector hull of J 1 /) is 
naturally identified with TM: 

= TM. 

Now we consider the bundle /?2,i : J 2 P ~^ J 1 / - We have seen in section 5.1 that it is an affine 
bundle modelled on the vector bundle Vpifl, which is a subbundle of the tangent bundle T(J 1 / o), 
and that J 2 p can be naturally embedded into T(J x p). Our aim is to identify the vector hull J 2 p 
with a suitable subbundle of T(J 1 / o). 
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The Cartan distribution Cpifi is a subbundle of T(J 1 p) which includes Vpi ; o, and it is 
straightforward to see that the sequence of vector bundle morphisms 



>■ Vpifl c >■ C Plfi — ^ J 1 p x R >■ 

is exact, where w denotes denotes the contraction with the 1-form dt, restricted to vectors 
in Cpifi. It is easy to see that is equal to J 2 p as a submanifold of T(J 1 p), so we can 

conclude that J 2 p is naturally identified with Cpip: 

J 2 p = Cpifl. 
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